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Fluctuation-dissipation theorems for viscoelastic fluids: 
consistency of extended thermodynamic theories 
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Fluctuation-Dissipation Relations (FDR) for a Maxwell fluid are computed via the GENERIC 
formalism. This formalism is determined by four building blocks, two "potentials" (total energy and 
entropy) and two "matrices" which determine the dynamics, but the understanding of fluctuations in 
a given non-equilibrium system is reduced to determining a single friction matrix. The FDR exhibits 
interesting features, arising from the type of entropy used in the formalism. We show explicitly this 
dependence FDR-Entropy, also we show that GENERIC renders results consistent with irreversible 
linear thermodynamics (TIL) only if it is used the corresponding entropy. An inconsistent result is 
provided by GENERIC when it is used a modified entropy . 
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I. INTRODUCTION 

Much effort has been expended within the theory of 
stochastic processes to obtain the mesoscopic basis of 
macroscopic theories. The problem of the description 
of equilibrium and non-equilibrium fluctuations in 
macroscopic systems is one of the central aspects of both 
thermodynamic and microscopic theories of irreversible 
processes 0. From the macroscopic viewpoint, the 
non-equilibrium case has been investigated within the 
framework of several theories among which we men- 
tion in particular extended irreversible thermodynamic 
theory (EIT) Q . The starting point of extended thermo- 
dynamics is the generalization of the Gibbs relation for 
the non-equilibrium entropy, which is used to determine 
the second moments of the physical fields under the 
assumption that the probability of the fluctuations is 
given by the Einstein relation |3|. Another, general 
formalism for studying the dynamics of non-equilibrium 
systems is called GENERIC (general equation for the 
non-equilibrium reversible irreversible coupling) 4] . This 
framework has been developed to formulate dynamic 
equations for out-of-equilibrium systems, by empirical 
arguments and by identifying the common structure 
of successful models in non-equilibrium thermodynam- 
ics and it has been validated by projection operator 
techniques [f|. In particular, Generic has proved to be 
a strong framework for dealing with the modeling of 
rheological properties of complex fluids • 
Recently, in an attempt to test the consistency between 
alternative theories of irreversible processes, the equi- 
librium Rayleigh-Brillouin spectrum for a Maxwell fluid 
has been computed. This was made by using fluctuating 
hydrodynamics, via fluctuation-dissipation relations 
(FDR), and by the procedure used by Mountain || 



which is a direct calculation technique. One would a 
priori have expected that since the Maxwell model has 
been often quoted as the prototype system for an EIT 
treatment, the equivalence of the two routes would have 
occurred when performing both types of calculations. 
However, a discrepancy arose whose origin could be 
traced to the fluctuation-dissipation relation for the 
(stochastic) stress tensor. If the prescription made by 
Landau for the frequency-dependent dissipative coeffi- 
cients case is used both routes coincide. In contrast, if 
the FDR from EIT is used the two routes diverge. 
In this work, we undertake again the testing of consis- 
tency between alternative thermodynamic theories. We 
describe fluctuations in a viscoelastic fluid by using three 
different thermodynamic frameworks, namely, linear 
irreversible thermodynamics (LIT), and two versions of 
extended irreversible thermodynamics. The first one is 
extended with the stress tensor, and the second one with 
the conformation tensor of the kinetic theory of poly- 
meric fluids [4|. Our purpose in this paper is to discuss 
explicitly, via GENERIC, why the LIT framework must 
to be used for calculating a FDR in rheological models 
for complex fluids instead of the extended versions of 
irreversible thermodynamics. An unavoidable conclusion 
is arising from this analysis. The extended versions of 
fluctuating irreversible thermodynamics must be revised. 
The paper is organized as follows. Section II is devoted 
to describe the main characteristics of GENERIC to 
make this article self-contained. In Section III we 
illustrate how the FDR are obtained for a known case, 
namely, a Newtonian fluid which is described by LIT. 
We analyze in Section IV the case of a viscoelastic fluid 
within LIT and EIT by extending the thermodynamic 
variable space with the dissipative fluxes and with the 
conformation tensor Q. We obtain the corresponding 
FDR for each case. Finally, in Section V we make 
some discussion on the previous issues and add some 
conclusions. 
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II. FLUCTUATION-DISSIPATION RELATIONS 
AND GENERIC 

The fundamental time-evolution equation of 
GENERIC, which was developed by considering 
the compatibility of two levels of description and by 
studying a large number of specific examples 0, can be 
written in the form 



dx SE(x) SS(x) 

-dt= L{x) — + M{x) — - 



(1) 



where x represents the set of independent state variables 
required for a complete description of the underlying non- 
equilibrium system, the real-valued functionals E and S 
are the total energy and entropy expressed in terms of 
the state variables x, and L and M are the Poisson and 
friction matrices (or linear operators). The two contribu- 
tions to the time-evolution of x generated by the energy 
E and the entropy S in Eq. ^ are called the reversible 
and irreversible contributions to GENERIC, respectively. 
Since x typically contains position-dependent fields, such 
as mass, momentum and energy densities, the state vari- 
ables are usually labeled by continuous labels in addition 
to discrete ones. A matrix multiplication, or the appli- 
cation of a linear operator, hence implies not only sum- 
mations over discrete labels but also integrations over 
continuous labels, and S/5x typically implies functional 
rather than only partial derivatives. 
In the GENERIC framework, Eq. J3> is supplemented 
by the complementary degeneracy requirements 



L(x)—^- = 0, 



8x 



SE(x) 
5x 



0. 



(2) 



(3) 



The requirement J5J) expresses the reversible nature of the 
L-contribution to the dynamics: the functional form of 
the entropy is such that it cannot be affected by the op- 
erator generating the reversible dynamics. The require- 
ment (01 expresses the conservation of the total energy 
by the M-contribution to the dynamics. Furthermore, it 
is required that the matrix L is antisymmetric, whereas 
M is Onsager-Casimir symmetric and semi-positive def- 
inite. Both the complementary degeneracy requirements 
(|2~|). J2Jl and the symmetry properties are extremely im- 
portant for formulating proper L and M-matrices when 
modeling concrete non-equilibrium problems 0. Finally, 
the Poisson bracket associated with the antisymmetric 
matrix L, 



{A,B} 



5 A T , JB 
— , L(x) — 
ox ox 



(4) 



with 



{A,B} = -{B,A} (5) 
is assumed to satisfy the Jacobi identity, 

{{A, B},C} + {{B, C},A} + {{C, A},B} = 0, (6) 

for arbitrary functionals A, B, and C. This identity, which 
expresses the time-structure invariance of the reversible 



dynamics, is postulated as another important general 
rqperty required by non-equilibrium thermodynamics 



The stochastic dynamics into GENERIC is determined 
by the stochastic differential equation obtained by adding 
noise (and the divergence of M) to Eq. .3, 

dx = L^ + M^- + k B ^dt + B{x)dW u (7) 
ox ox ox 



where B is a solution of the equation 



BB T = 2k R M. 



(8) 



and Wt is a multi-component Wiener process, that is, a 
Gaussian process with first and second moments given 
by 

The expression (jHJ for B may be regarded as the 
fluctuation- dissipation relation. 

We now show how the FDR may be obtained in this 
framework. Eq. Q is written by components as follows 

dx?(t) = A?{x)dt + B™™' (x)dW™' (t), (10) 

where the definition of A™(x) is evident and the incre- 
ments dW satisfy the property 



dW™(t)dWp (t) 







if t ^ t' 



(11) 



The Kronecker's delta S rr i is equal to unity if r = r' 
and vanishes otherwise. In Eq. I|1U|) a white noise term 
defined as 



rT(r,t) 



= dW™(t) 
~ dt 



(12) 



may be introduced. This term has the following correla- 
tion function 



■rMwy.o)- ^"^ ) 

I $mm' $rr' if t t , 

jO if t ^ t'. 



(13) 



This expression may be rewritten as 

(»T(r, t)rT' (r', t')) = 5 mm ,8{r - r')S(t - t'). (14) 

which gives the same results as Eq. (|13|l . We define the 
stochastic term in Eq. IjlOjl as 

f?{t)dt = ^BZr'{x(t))dWF\r,t). (15) 

r' 

Calculating its correlation we obtain 
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dt / dt"f?(t)fp, (*") 



/ / M*))^" 1 ' t) x ^ B™'' r w"' {x{t'))dW m "' (r'", t') 

dtJ^B™r \x{t))B^ \x(t))\ 
2k B j dt j Dxf(x,t)M™™"(a) = 2k B [ dt(M™? 



(16) 



r 



where we have made use of Eq. (Illf) . f{x,t) is the prob- 
ability density that the system takes the values a at time 
t Moreover, note that 

/™(t) /#"(«")) = 2fc B *(t - t") (AO?") 4 (17) 

since if we integrate on i, t" we recover Eq. 1)160 and 
the integration has arbitrary limits. Expression (|17|) is 
another way to write the FDR. We will use this expres- 
sion in the following Sections. As it has been mentioned 
the form of matrix M depends on the form of the en- 
tropy out from equilibrium. Therefore it determines the 
specific form adopted by expression (|1 7|) . 



III. THE VISCOUS FLUID IN LIT 

As an illustration of the stochastic dynamics and 
, we consider an example of fluctuating hydrodynam- 
ics. The results of the previous Section are illustrated 
with the case of a viscous Newtonian fluid within LIT. As 
usual the physical conserved densities are chosen as inde- 
pendent variables to describe the thermodynamic state of 
the fluid. Mass density is denoted by p, momentum den- 
sity by u, and internal energy density by e. The balance 
equations of mass, momentum and energy are written as 



d 8 



du d d d 

dt dr dr dr 



(18) 



(19) 
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(20) 



de d d d 

dt dr dr dr dr 

where p is pressure, the expression for the stress tensor 
is that of Newton and it is written in terms of the veloc- 
ity gradient, shear viscosity 77, and volume viscosity £ as 
follows 



d 

— i 

dr 



d 

— i 

dr 



c 



d_ 
dr 



vl (21) 



Fourier's heat conduction is assumed and written in 
terms of the temperature gradient and thermal conduc- 
tivity K 



dr 

dr 



(22) 



Knowing the temperature T = T(p, e) and pressure 
p = {p,e), Eqs. (|18I22I) form a closed set of time evo- 
lution equations. The irreversible contributions to the 
dynamics are those involving the stress tensor *r* and the 
heat flux q appearing in Eqs. (fTT)|l y (J20J) ; the remain- 
ing terms are the reversible contributions to the dynam- 
ics. The dissipative coefficients may depend on position. 
Though the matrix of the reversible part of the dynamics 
is not necessary for the study of fluctuations we include 
a brief summary to give the complete GENERIC'S equa- 
tion. The total energy needs to be expressed as a function 
of the set of variables of the system. This set is denoted 
by the vector x — (p, u, e). The vector SE/Sx results in a 
vector with five entries constituted by functional deriva- 
tives of the energy with respect to x. The total energy 
is obtained by integrating the addition of the kinetic en- 
ergy and internal energy densities over the volume of the 
system 



E = 



l u(r)' 
2 p 



e(r) 



d 3 ? 



(23) 



The functional derivative of this energy with respect to 
x, is given by 



SE 
Sx 



Sp(r) 

J^r) I E(p,u,e) 



Se(r 




(24) 



The matrix L may be constructed in two ways. Both 
have been explored in detail in Q and ja| • We only men- 
tion here that the last authors have based their method 
on variational techniques. To determine the dissipation 
operator, the functional of entropy of LIT is considered 



(25) 



S= / s(p(r), £ (r))^r 
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in such a way that the functional derivative with respect 
to the hydrodynamics fields is given by 



SS 




— - . « . (26) /i(r) dsfc.e) 

fell) ^y= ^- (28) 



where the local temperature and the chemical potential 
are defined as 



ds(p, e) 
de 



The matrix M, which reproduces the irreversible terms 
1/T(r) (27) in the dynamic equations, then takes the form 

I 



M s (r, r') 











' d d 
\ dr' dr 



d d 
dr dr' 



d_ 



& ■ V T 7 S + £rkT(tv[j])S 



■rfTtf + &kT(tr[y])8 



(29) 



where 



7 = 



d < ^ 



2 3 



Integrating by parts over r', we get the matrix M which 
depends only on r 



(30) 



M s (r) 



(0 



V> 







_ T? r 7 &-fcT(tr[ 7 ])& 



9 9 a 

— l\i _i TT — 

ar ar 



5 7 ? T 7 : 7 



(31) 



9 ■ ^T 2 # + |feT(tr[ 7 ])V 



,7r 



The elements of this matrix are differential operators act- 
ing on the terms on its right including the vector 5S/6x. 
the matrix M s (r,r') satisfies the required degeneracy 
condition 



M ■ 



5E 
Sx 



f M s (r,r')- 1 




v(r) 1 dV = (32) 







The physical interpretation of this condition has been 
given in Section II. We now calculate the FDR for the 
stress tensor by using the last expression for the dissipa- 
tive operator. A similar procedure may be followed for 
the heat flux. We omit it in this work for simplicity. The 
FDR is given by Eq. (|17|) . As the stochastic term asso- 
ciated with the momentum equation is introduced as a 
divergence, which we denote by V*, we write firstly the 



FDR as follows 



ir-f (*) 4r<> (*': 



di 



= 2k B 5(t~t')^C) 

(33) 



where the term M^, is determined by the corresponding 
elements of the dissipation operator matrix Eq. I|29|l . We 
have explicitly 



M -=(4^ + *"^r r)T(r,M - r ' ) 

(34) 
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which may be rewritten as 
rr dr n dr 



Mt =— — {{b av b^ + mf>») r) (r) T (r) 



+ 2 Qc W - \r, (r)) T (r) S aP S^}6 (r - r') 



(35) 



Therefore the FDR takes the form 



(af (t) o£" (*')) = k B {(5 av 5^ + 5 a »5^) (2r, (r) T (r)) 



(36) 



+ 2^((r)-- )) (r)jT(r))x 
5 a p5^}5 (r-r') 5 (t-t') 



This result coincides with that of Espahol 0] obtained 
from a projection operator technique. It corresponds to 
a viscous fluid out of equilibrium and the dissipative co- 
efficients C arL d rj are dependent on the position. When 
the dissipative coefficients are constant and the system is 
in equilibrium the mean value of the temperature is the 
equilibrium temperature Tq and the FDR becomes 



{af (t) o»r (f)> = k B {{S au 6^ + S a ^ v ) 2 V T 



S a/3 S^}S(r~r')S(t-t') 



(37) 



reproducing the Landau-Lifshitz's FDT for the stress ten- 
sor [13. 

IV. THE MAXWELL FLUID 

In this Section we deal with a viscoelastic fluid and 
study the fluctuating properties of this kind of fluids. 
One way to study viscoelastic fluids from a thermody- 
namic point of view is by assuming that viscosity and 
elasticity can be considered separately and that the net 
effect is additive. Under this assumption one arrives to 
one of the simplest models for viscoelasticity, namely, the 
Maxwell model. A complete derivation of the Maxwell 
model in the macroscopic level may be seen for example 
in . One then may tackle the problem of the fluctu- 
ations by inserting the Maxwell model in different ther- 
modynamic frameworks. We concentrate here on three of 
them. On the one hand, LIT in which the state variables 
of the system are the usual conserved densities of mass, 
momentum and energy. In this formalism the entropy is 
considered a function of the conserved densities solely as 
mentioned in the previous section. On the other hand, we 
consider extended theories. Different ways of extending 
the state variables space in order to deal with viscoelas- 
tic fluids exist giving rise to different levels of descrip- 
tion. A macroscopic level of extended thermodynamics 
is obtained when one extends the variable space with the 



dissipative fluxes of the system, the heat flux and the 
stress tensor in this case. We get one more microscopic 
level when the same space is extended with structural 
variables giving an explicit treatment to some internal 
degrees freedom of the constituents of the fluid. In this 
level of description one must assure that the dynamic 
equations of the more microscopic state variables repro- 
duce the Maxwell equation for the stress tensor. In both 
of these extended versions of irreversible thermodynam- 
ics the entropy is considered a functional of the enlarged 
variables space and this introduces a rather controver- 
sial issue in the problem. Thus the system under study 
is a viscoelastic fluid modelled as a Maxwell fluid when 
the relaxation times of the fast properties of the system 
(dissipative fluxes, structural variables) are of the order 
of the observation time. One would expect that different 
thermodynamic frameworks lead to the same physical de- 
scription of the system. Surprisingly, the results diverge. 
We discuss, in some detail (the final comments in Section 
V), how to discern the validity of these diverse results by 
considering the light dispersion properties of the fluid as 
obtained from the different thermodynamic frameworks 
investigated. 

We first examine the LIT case. This problem must be 
treated in Fourier space so we firstly write the GENERIC 
equation in that space. Consider then the next notation 
for the time Fourier transform 



/M 



i 



J- 



f (t) exp (iwt) dt 



Then the GENERIC equation becomes 



-in7a; (vj) 



~ SE 
L * — 
ox 



~ SS 
M * — 
ox 



(38) 



where the symbol * denotes the convolution of the in- 
volved Fourier transformed functions. It is defined as 
follows 



f*9 = 



1 



f(t')g(w-t')dt' 



(39) 



Given the dynamic equations of the system the problem is 
then to find matrices L and M in the Fourier space which 
substituted in Eqs. (|3*5|) lead to the dynamic equations 
of the system. We now write the linearized equations of 
the Maxwell fluid in the Fourier space 



d „ 

-impx = -Pott ■ vi 
or 



-impovi 



dp 3 p; 
dr dr 1 



-iwCTi 



-p * 



dr 



d_ 
Or 



dTi 
dr 



(40) 



(41) 



(42) 
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These are the dynamic equations in the variable thermo- 
dynamic space of LIT. They must be supplemented with 
the time evolution equation of the Maxwell-Cattaneo 
kind for the stress tensor 



( I - iz^T, ) Vj = -2// ( ^-Vl j - \ //, 



V x 



(43) 



Eqs. H4UI43[I have been linearized around the equilib- 
rium state. So we have used the following approximated 
expressions for the dynamic variables 

P = Po + Pi 



By substituting Eq. I|43|) in Eq. (|41|l the momentum 
density balance equation takes the form 



dp d 
dr dr 



dr 



vi + 



where the dispersive coefhcients are given by 



D I \ (1 + i^Tr) 77 



1 + m 2 T? 



(44) 



v = vi 



T = T +T 1 



T = T 1 



K (w) = 



(1 



(Ha - H) 
v 2 3/ 



We now concentrate on the dissipative part of equations 
gnjl, (ISI) and igU and propose the next matrix M 



M = 



(& £r+^l§z-£r)v D (^)T 8 + 2££rK (m)T 5 2& ■ rj D (w) T 7 8 + (m)T Tr j 5 



\0 ( w )T Tr 1 5 



r 



(45) 



It is a direct task to verify that this matrix when applied The resulting FDT is 
on the first order divergence of the LIT entropy 




reproduces the dissipative part of Eqs. ^4Up. (|44|l and 
(|42|1 . Now we use the expression Eq. I|33|) to find the 
fluctuation-dissipation theorem for the stochastic part of 
Eq. (|44|l . which is written as a spatial divergence as in 
the Newtonian fluid. The central block of Eq. (|45|l de- 
termines the form of the fluctuation-dissipation theorem. 
The form of the FDT in the Fourier space is 



'_d_ 

>dr a 



^M*^(-')) 



(af H * 5£" (a/)) = k B \2 (8 a »8^ + 8 a "8 0v ) x 



t] D (zu) T + 25 al3 S^K (ro) j x 

8 (r - r')S {w + tu') (46) 



which corresponds to the prescription for the FDT made 
by Landau |12| in the case of dispersive transport coeffi- 
cients. 

We now examine the case of the extension of the variable 
space with the dissipative fluxes, specifically, with the 
traceless stress tensor and its trace in such a way that 

x = ^p,v,e,T,T^j , we define the state of equilibrium as 

x o = {Po, 0, £ , 0, 0, 0). The balance equations, linearized 
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for simplicity, are given by 



d d . 



(47) 



de d . . d d d 

«7 = --^—( £ oV a ) -p-^—v a + -^—-^—KTi 

at or a or a or a ar a 

and the relaxation equations for the extended part of the 
variable space are 



T 2 - 



dt 



-2t]X^ pa — v p (48) 



dr 



The entropy of the fluid is assumed to be quadratic in 
the fluxes 



S{p,s,T a0 ,T) = S eq - — r — 4^T Tq/3 : T ! 



,a/3 



Calculating the functional derivative of the entropy with 
respect to the traceless stress tensor and its trace we get 



5S 



St 



a/3 



VT 2 o 



VT2 o 1 



T 



ss 
5^ 



VT 1 



(49) in such a way that matrix M takes the form 
I 



(51) 



(52) 



/0 




-^O-rSSuufl 

p ACv or a or fi ^ u ° 



4:1] 



M(r, r') 











v'T a (, tA- 1 



V : ) O 



We note that the two last entries of the vector 5E/8x 
vanish since the total energy does not depend on the 
traceless stress tensor and its trace. Therefore, the aug- 
mented matrix M, Eq. (|53[) . trivially satisfies the degen- 
eracy condition on the gradient of the energy. Eq. I|17() 
does not lead to the right result for the FDR. See for 
example Q, where stochastic process theory was used 
to obtain the FDR for the stress tensor of a viscoelastic 
fluid modelled as a Maxwell fluid. One may also prove 
that the quadratic form of the entropy, Eq. (|5U|) does not 
satisfy the degeneracy condition, Eq. J5| with L given in 
for the viscous fluid, making this description thermo- 
dynamically non-consistent. 

We now describe the viscoelastic fluid within the ex- 
tended version given in the original GENERIC'S papers 
P| . The variable space is enlarged with the conformation 
tensor c given by 



c(r) 



1 



QQV'(r,Q)d 3 Q 



(54) 



where ip(t, Q) is the configurational distribution func- 
tion, Q is interpreted as a dumbbell configuration vector, 
and rj p is the polymer concentration considered here as 



constant [4|. By considering that the complex character 
of the viscoelastic fluid has only entropic consequences 
the total energy of the fluid is written as follows 



E 



lu(r) 



+ e(r) 



d 6 r 



2 p 

While the entropy functional becomes 

S = j s(p{r),E(r))d 3 r + S p 



(55) 



This expression includes the polymeric contribution to 
the entropy in the term S p . This term is given by |4| 

S p = -VpkB J {tr [l-cc(r)] + In [dct cc(r)]} d 3 r (56) 

The osmotic pressure tensor is chosen such that the en- 
tropy gradient belongs to the null space of the Poisson 
operator 



n = T[2c-^ + 5 p l 
dc 



(57) 
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Substituting Eq. in Eq. ({5T|) it is obtained 

S = r] p k B T{l - cc) (58) 
I 

M = 

In this way the form of matrix M, together with the 
matrix L (which is very similar to Eq. (54) of ||) repro- 
duces the Maxwell equation for the time evolution of the 
osmotic pressure tensor 

( 1 + A "^)v^T = - « (60) 
I 



The matrix M, was proposed by Ottinger and Grmela in 
0, as follows 



(59) 



I 

where D/Dt is the material derivative. However this 
form of M is not suitable to be applied to obtain the 
FD relation since the lower right term does not have the 
proper tensor rank. A simple possible modification of M 
is 



/0 \ 

~VTi£ hTl : 7 - & • KT^JL 

V° ^CT C / 



M = 








Or / 







(61) 



r) p k B c\ H 



-lc 



where we have changed c- by lc : . In this way, the degen- 
eration condition Eq. J3J is preserved and the Maxwell 
equation (|6T7|) is again obtained. 

However we make the reader note the fact that this last 
form of M leads to a FDR which is in discrepancy with 
that expected for the osmotic component of the stress 
tensor in the case of a viscoelastic fluid 

0.0- 

Never- 
theless it must be mentioned that the obtained FDR is of 
the (^-correlated type. This fact will be used in our final 
discussion. 



V. DISCUSSION 

GENERIC 4] is a formalism which describes the dy- 
namics of complex fluids in terms of two generators, the 
total energy and the entropy of the system. The thermo- 
dynamic consistency of the time evolution equations of 
the variable space of the system is assured through the 
degeneracy conditions, Eqs.J2J) and ©, on the matri- 
ces containing the conservative and dissipative dynamics 
of the system. We remind the reader that the physical 



meaning of the two mentioned conditions is the second 
and first laws of thermodynamics respectively. The ther- 
modynamic consistency of any constitutive model is then 
assured if the two conditions are satisfied by the con- 
stituents of the GENERIC form of the dynamics equa- 
tions of the system. We also remark the key role played 
by entropy in determining the explicit expressions for the 
fluctuation-dissipation relations in the fluctuating hydro- 
dynamics obtained from the GENERIC formulation of 
the dynamics. 

We have used the fluctuating formulation of GENERIC 
to find the FDR for a Maxwell fluid in each one of the 
three thermodynamic frameworks considered. On the 
one hand in Section IV we have already given the rea- 
sons why the EIT version with the variable space en- 
larged with the dissipative fluxes does not lead to valid 
results. We remind the reader that the expression for 
the entropy which leads to the dynamic equations for the 
Maxwell fluid does not satisfy the degeneracy condition 
Eq. (2). In their papers Grmela and Ottinger also noted 
that the quadratic form of the entropy for the Maxwell 
fluid does not satisfy the condition (J2J. Specifically, if 
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one would assumes the quadratic form for the polymeric 
entropy 

S P = —Vpka J [l-cc(r)] : [l-cc(r)] d 3 r, (62) 

then the gradient of entropy (S = S eq + S p ) would not 
be in the null space of the Poisson operator. Thus, the 
quadratic form of entropy is not a suitable one from the 
point of view of its thermodynamic consistency. 
We now concentrate on LIT and EIT version enlarged 
with the conformation tensor given in Section IV. In or- 
der to discern which of these results is valid we resort to 
a previous work by Vazquez and Lopez de Haro §. In 
that paper the equilibrium Rayleigh-Brillouin spectrum 
for the Maxwell fluid was computed. This was done by 
using a procedure by Mountain and by fluctuating hy- 
drodynamics based on two types of FDR. The first of 
them was of the (5-correlated type and the second one of 
the exponentially correlated type which is the one pre- 
scribed by Landau and Lifshitz for the case where a dis- 
persive dissipation coefficient is present. In the first case 
the Maxwell fluid was described in the EIT framework 
and in the second one the fluid was described within 
LIT. Mountain and fluctuating hydrodynamics calcula- 
tions rendered the same result only if the Maxwell fluid 
was described within LIT and the corresponding expo- 
nentially correlated FDR was used. From this fact one 
should conclude that a 5-correlated FDR does not lead 
to the correct description of the Rayleigh-Brillouin spec- 
trum. This is the case of the FDR in the extended version 
of the viscoelastic fluid GENERIC theory. At this point 



we might conclude a number of facts in spite the situ- 
ation is rather unsatisfactory. Firstly, we remark that 
the GENERIC scheme allows an effective definition of 
the entropy of the system by respecting two fundamen- 
tal physical principles, namely, the first and second laws 
of thermodynamics. The EIT extended with the dissi- 
pative fluxes leading to the quadratic form for entropy 
is disqualified by its thermodynamic inconsistency. The 
EIT extended with the conformation tensor will not ren- 
der the expected Rayleigh-Brillouin spectrum since the 
Maxwell fluid's FDR within the formalism takes the form 
of a (5-correlated FDR. Secondly, LIT is consistent from 
two points of view. On the one hand, it satisfies the 
prescribed degeneration conditions of GENERIC and on 
the other, it describes in the correct way the fluctuating 
properties of the Maxwell fluid leading to the expected 
Rayleigh-Brillouin spectrum. Thirdly, thus one is forced 
to conclude that the extended versions of fluctuating irre- 
versible thermodynamics should be revised. We end this 
paper with the following conclusion. LIT is a proper ther- 
modynamic framework to deal with even situations out 
from equilibrium for which the extended theories have 
claimed to require a set of variables in addition to the 
usual set of conserved densities of LIT to characterize 
the thermodynamic state. 
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